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Abstract 

Raman scattering plays a key role in umaveling the quantum dynamics of 
graphene, perhaps the most promising material of recent times. It is crucial to cor¬ 
rectly interpret the meaning of the spectra. It is therefore very surprising tiiat the 
widely accepted understanding of Raman scattering, i.e. Kramers-Heisenberg-Dirac 
theory, has never been applied to graphene. Doing so here, a remarkable mecha¬ 
nism we term ^Transition sliding’’ is uncovered, explaining the luicommon brightness 
of overtones in grapheue. Grapheue’s dispersive and fixed Raman bands, missing 
bajicls, defect density and laser frequency dependence of band intensities, widths of 
overtone bands, Stokes, anti-Stokes anomalies, and other known properties emerge 
simply and directly. 


1 Introduction 

The imique properties of graphene and related systems have propelled it to a high level 
of interest for more than a decade. Raman scattering is perhaps the key window on 
graphene’s quantum properties, yet the crucial aspects of the spectriun of graphene, cai- 
bon nanotubes, and graphite Raman spectra have been a subject of much controversy, in 
some cases for decades. 

The delay in applying Kramers-Heisenberg-Dirac theory (KHD) to graphene is at¬ 
tributable to the rise of a non-traditional, 4th order perturbative Raman model (KHD 
is 2nd order) called ‘'double resonance” (DR), which appeared 15 or so years ago[I], 
rapidly gaining wide acceptance exclusively in the conjugated carbon Raman commu¬ 
nity. Far from being a variant of KHD tuned to graphene, DR is starkly incompatible 
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with KHD. DR. exclusively relies on post-pliotoabsorption inelastic electronic scattering 
to create phonons, requiring introduction of two additional orders of electron-]>honon per¬ 
turbation theory^ making DR overall fourth order. These make no appear*ance in KHD, 
which cleanly produces phonons free of Pauli blocking by two other mechanisms, overall 
in second order perturbation theory. 

Search engines quickly reveal that none of the 3000+ published contributions devel¬ 
oping DR. or using it to interpret experiments has mentioned the much older and more 
established KHD Raman formalism, in the last 15 years. No rationale was ever given as 
to why a radically different, non KHD formalism was needed. On the face of it, even with¬ 
out regard to the tlieory given in this paper, these facts leave the DR approach seriously 
unexamined and exposed. We show here that imfortunately DR. was in fact a wrong turn 
on the way to understanding graphene quantum dynamics and Raman spectroscopy. 

When KHD is a]>plied to graj^hene, phonons are produced exclusively and instantly 
at the moment of photoabsorption or photoemission by the long neglected (in condensed 
matter physics) nucleai* coordinate dependence of the transition moment. Graphene’s 
Raman spectroscopy falls quickly and naturally into place, providing new insights and 
predictions along the way. 

The present graphene study has a strong precedent in another conjugated conducting 
polymer. In a recent study [2], the long-standing mysteries of the dispersive polyacety¬ 
lene Raman spectrum yielded to KHD theory, extended to include 1-D crystal structure, 
defects, and electron and phonon dispersion relations. All these arise from within Born- 
0]>peuheiiner theory and KHD Raman scattering theory. The coordinate dependence 
of the transition moment plays a key role, since geometry does not change when the 
electron-hole pairs are created in the extended conjugated system. 

These principles are the starting point for surprising new insights that graphene has 
in store. For example, '‘transition sliding’, not possible in a DR. picture, is responsible 
for the brightness of overtone bands, and explains their impurity and doping dependence. 
Transition sliding does not apply to polyacetylene because there are no Dirac cones; indeed 
polyacetylene overtones ai*e correspondingly weak. 

2 Virtual processes in KHD versus DR 

DR makes heavy use of virtual states that are foreign to ABO. ABO is certainly not exact 
(see next section), so for example its predicted electronic band structure is not exact. Still 
it is siuprising to see virtual states freely used in DR. that do not exist on the electronic 
bands, or virtual states that violate momentum conservation, or even violate the Pauli 
exclusion principle [ 3 . 

If one accepts ABO as at least close to the truth, no such virtual states can arise. 
All the terms in the KHD siun are inomeutum conserving, Pauli principle obeying states 
lying exactly on Dirac cones or other parts of the electronic dispersion surfaces. KHD is 


2 


based on second order time dependent light-matter perturbation theory applied to Born- 
0]>penheimer states. Only the energy, conjugate to time, can become virtual, as familiar in 
off-resonant Raman scattering. Even “temporary" violation of momentum conservation 
is forbidden in KHD, since momentum violating terms in the sum in equation give 
vanishing matrix elements, whether they are off-resonance or not. 

3 Does the Born-Oppenheimer Approximation “Break 
Down” for Graphene? 

KHD relies lieavily on tlie Born-Oppenheimer approximation, as does much of condensed 
matter physics. Recently, the “headlines^’ in a highly cited paper left the impression that 
the ABO fails particularly in graphene. The article was excellent, article with the provoca¬ 
tive title "'Breakdown of the adiabatic Born-Oppenheimer approximation in graphene^' J3]. 
published in Nature Materials with over 800 citations, we find the statements ABO 
has proved effective for the accurate determination of chemical reactions, molecular dy¬ 
namics and phonon frequencies in a wide range of metallic systems. Here, we show that 
ABO fails in graphene." and later “Quite remarkably, the ABO fails in graphene” [3], see 
also [ 3 . Such statements read like an announcement of catastrophic failure of ABO in 
gi'aphene. Several times in the course of explaining the present work, one of the authors 
has been asked, ‘"but didn’t you know that the Born-Oppenheimer approximation fails in 
graphene?” 

Of coiurse the Born-0]>penheiiner approximation fail in graphene, and in every other 
system it has been applied to in the last 90 years, in specific regimes and situations. It 
is, after all. an approximation. Reference J3] claims nothing so drastic as a catastrophic 
break<lown once one reads deeper. Instead it makes important points about a stiffening of 
the G mode with electron density near the Dirac point involving the Kohn anomaly, missed 
by ABO. By their nature Kohn anomalies involve rapid changes of electronic structure 
with small nuclear configuration changes, a warning sign of potential ABO breakdown. 

In this paper we do not attempt to correct for the Kohn anomalies and their effect on 
phonon modes and mode frequencies. Electron-hole pairs created in Raman studies are 
normally well away from the Dirac point and orbitals affected by Kohn anomalies. 

The Born-Oppenheimer approximation, proposed in 1927, is the basis of crystal struc- 
tiure. electronic band structure, and phonons and their band struct me. It is quite arguably 
the main pillar of condensed matter physics. Surely it needs monitoring and corrections, 
but without it we have a structureless pea soup of electrons and nuclei, rather than crystals 
and solids. 

The ABO succeeds in the large, that is the point. Recall that within ABO, if the 
nuclei return to a prior configuration, the electrons do also. This is strictly incompatible 
with inelastic electron-])honon scattering having taken place in the meantime. Inelas¬ 
tic electron-phonon scattering is nonetheless a valid non-adiabatic correction to ABO. 
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However, electron-phonon scattering would appear to Pauli block the affected electrons^ 
immediately rendering them silent in Raman scattering. 

We have been given a late opportunity to use KHD theory the first time in graphene. 
ABO is vastly more established than DR, and astonishingly accm’ate most of the time. It 
much simpler to execute, and lower order in perturbation theory. We hope the community 
is curious as to how it performs. 


4 Kramers-Heisenberg-Dirac Theory 

We provide a review of KHD here: facilitating the developments in succeeding sections. 
Just before the dawn of quantimi mechanics, in 1925 Kramers and Heisenberg published a 
correspondence principle account of Raman scattering which Dirac translated into 

quantum form in 1927 jHj. The Kramers-Heisenberg-Dirac Theory of Raman scattering 
has been used ever since to explain more than half a million Raman spectra in a very 
wide range of systems, including very large conjugated hydrocarbons. There is no reason 
that removing hydrogens from cai'bon materials should cause KHD to catastrophically 
fail and require replacement by a theory based on different physics. 

The KHD formula for the total Raman cross section E reads, for incident frequency 
oj[ and polarization p, scattered frequency and polarization cr, between initial Born- 
0]>penheimer state \i) and final Born-Oppenheimer state |/) via intermediate Born- 
Oppenheimer states |n) reads 

(/|Dt-^|n} (n|D^|i} (/|Dt>^|n} (n|D^ 

Ei — En + tuJi — iTn Ei + En + tiojj + 

( 1 ) 

where is the damping factor for the excited state , accoiuiting for events and degrees of 
freedom not explicitly represented. The transition moment operator D controls the first 
order perturbative matter-radiation coupling. Usually the second, non-resonant term 
inside the square in equation is neglected, for simplicity. 

Making this more explicit, suppose with phonon coordinates £ and electron coordinates 
r, we write 



E?:!. = 


9c^ 
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l<>l 




I/) = \<?-{Cr)) = m-.r))\Xr.m^ 

N) - K°(Cr))^|^(Cr))|Xm,(0)>aiKi 

r)) is the approximation to the Born-Oppenheimer ir electron ground state (we su]>- 
press the other electrons) based on a Slater determinant of valence electron spin orbitals; 
l^qtqc(C*^)) electron-hole pair relative to the groimd state, with an electron in the 
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conduction band orbital |«Pqc(€:r)) and a hole in the valence band orbital | 0 q^,(^:r)). 
l^qtqc(C* ^))IXinvc(C)) ^ coiuplete intermediate state, including the phonon wavefunc- 
tion. q^,, and niyc niay range quite freely, with the following remarks: ( 1 ) almost all 
qi,qc. niic will give vanishing matrix elements with the dipole operator, due to momen¬ 
tum non-conservation (including the momentum of the phonons).(2) Some pairs qrqr with 
unmatched momentum q,. 7 ^ qr give non-vanishing matrix elements nonetheless because 
they are associated with newly created or destroyed phonons contained in |Xmvc( 0 ) rela¬ 
tive to the initial state, conserving momentum. Such states are the electron-hole-phonon 
triplets. (3) Momentum is conserved when elastic backscattering is required to re-align 
the hole and particle (because of the kick given to the conduction electron by newly 
formed or destroyed ]>honons) through the recoil of the whole sample, or laboratory, as is 
familial* in Mdssbauer scattering or elastic neutron scattering for example. This feature 
is embedded in equation since the electronic orbitals include the presence of impiuities, 
edges, etc., which are part and parcel of ABO. Even though the eigenfunction orbitals 
thus include the elastic backscattering as a boundary condition entirely within ABO, an 
electron promoted with a given momentum becomes a linear superposition of such states, 
and still takes time to backscatter. In the energy domain, being off resonance in effect 
shortens the time and nearly removes the effects of the backscattering. 


4.1 The transition moment operator 


We may write 


(n.|D'’|i) = (Xm..(C}|(K.qc(C;r)|D^|0(C;r))|xm{O). 


(3) 


as 

with 

- (0^.(tr)|D^(r)|0,,{^;r)),. (5) 

The matrix elements of the dipole operator D between two Born-Oppenheimer electronic 
states is the t7unsitio7i moment q^(€) connecting valence level qv and conduction band 
electronic levels qc. ^iqc,qv (0 written for light polarization p; the subscript r indicates 
that only the electron coordinates are integrated. Note that if it does not 

vanish, is explicitly a function of phonon coordinates 

In terms of the transition moments, the Raman amplitude reads (using only the res¬ 
onant term of the two) 


ail 
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The sum over the intermediate states witli energy denominators is of course the en¬ 
ergy Green function, i.e. the Fourier transform of the time Green function propagator. 
Although we do not write out the time dependent expression here (supplements, sec¬ 
tion 12.1 1 , we see that the valence state phonon wave function |Xm,(C)) in the 

conduction band multiplied and thus modified by the transition moment q(€).. includ- 
mg its phonon coordinate dependence, before any time propagation takes place. Moreover, 
that propagation does not further change the phonon pop^dations. 

The product of the transition moment and the phonon wavefunction is called 9qj,q,m, (0 
above. When (^qj,q,rn,(€) is re-expanded in all the phonon modes \Xm) conduction 

bands, some in will differ in the occupation luunbers compared to the imniodified initial 
lattice occupation (recall that because there is no geometry or force constant change, 
the phonons are the same in all valence and conduction bands). Every electron-hole 
pair production event carries some amplitude for no phonon change, plus si¬ 
multaneously a finite amplitude for phonon creation or destruction relative to 
the initial phonon wave function. 

The transition moment’s dependence on nuclear separation or equivalently phonon 
displacements is unquestionable, not only on direct physical grounds and explicit calcula¬ 
tions. but also because there would be no off-resonant Raman scattering without it. The 
derivative of the polarizability with nuclear coordinates, i.e. the Placzek polarizability, 
would vanish without coordinate dependance of the transition moment. If the transition 
moment is (unjustifyably) set to a constant independent of phonon coordinates, no Raman 
scattering occm’s at all in KHD theory for graphene. 

Due to dilution of the delocalized tt orbital amplitude over the infinite graphene sheet, 
it is almost obvious that the Born-Oppenheimer potential energy surface is unchanged 
after a single electron-hole pair excitation from the ground electronic state. No new forces 
on the nuclei arise upon a single electron-hole pair formation, and no geometry changes 
take place. Thus the independence of |Xmvc(€)) electron-hole state denoted by q^qr 

(the occupations indicated by vc may vary, but the phonon wavefunctions themselves do 
not depend on the electron-hole pair created). 

The intermediate, typically conduction band (assuming no doping) Born-Oppenheimer 
states IV^qJq^; *")) witli energy Ev,c,in range over resonant and non-resonant (very small 
or not so small denominators, respectively) states including all with non vanishing matrix 
elements in the numerator. Only those with Pauli and momentum matched electron- 
hole pairs and electron-hole-phonon triplets can be nonvanishing; it is not necessary that 
they be '‘resonant’*, i.e. minimize the denominator, in order to contribute importantly. 
Rather, a range of states are collectively important, the range depending on T. Although 
not relevant to graphene, even when none of the denominators reach resonance, there is 
still sufficient Raman amjditude to be quite visible as off resonance Raman scattering. 
The final Born-Oppenheimer state |\m/(0) with energy Emj conserves total system plus 
incident and emitted ]>hoton energy, and may differ from the initial state only by zero, 
one, two, ... phonons. The electron has filled the hole and the initial electronic state is 
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restored. 

All exception to this uotation is the filling an empty level created by doping, leaving 
the hole created by the photon unfilled. All the electronic states reside exactly on the 
electronic band surface, here on the Dirac cone neai* the Dirac point. 

It is good to keep in mind that the vast majority of states |0qyq,.(€: i‘))|Xnivc(€)) (vio¬ 
lating inoinentum conservation with respect to the initial state for example) in the KHD 
sum have vanishing matrix elements and do not lead to Raman (or even Rayleigh) scat¬ 
tering. Some terms (or really a small group of terms - see discussion in section about 
backscattering ) give non-vanishing matrix elements and if they include phonons, they 
are usually still incapable of Raman emission. An example is a D phonon produced in 
absorption in a clean sample, unable to emit a Raman photon due to the Pauli blocking 
effect of the electron recoil. Such states eventually relax and thermalize by a cascade of 
e-e and e-]>h inelastic processes. Transition (see section makes this worse in that even 
elastic backscattering by defects does not alleviate the Pauli blocking for the majority 
D phonons produced. These ‘‘wasted*' phonons, not appearing in the Raman spectrum, 
are always in the majority compai'ed to “successful” phonons. It is not the phonons 
themselves that fail, but rather their associated, Pauli-blocked conduction electrons. 

Even if no phonon is produced at the time of absorption, the chances of achieving 
Rayleigh emission or phonon production upon emission are small, given the withering 
removal of electrons from eligibility by inelastic e-e scattering, raising the near certain 
specter of Pauli blocking. This occurs on a femtosecond timescale (see pi fTHlfTTl [T^ fT^ ). 


4.2 Phonon adjusted electronic transitions 


A fraction of resonant conduction band eigenstates |0qvqc(€«r))|xmvc(O) one or 

more phonons relative to the valence band state. The energy of the Born-Oppenheimer 
eigenstates is a sum of phonon and electronic components, so the largest, resonant terms 
in the KHD sum equation necessarily arrive with the electronic transition energy cor¬ 
respondingly lowered if a phonon is activated in the conduction band eigenstate relative 
to the valence state or raised if a phonon is de-activated in the conduction band eigen¬ 
state relative to the valence phonon wave function. This key fact, missing in DR, alters 
electronic transition energies and phonon frequencies for dispersive bands. For phonon 
creation, the valence energy of the hole is raised, and the conduction band energy low¬ 
ered, keeping q the same in both, thus creating momentum matched (possible except 
needing elastic backscattering) holes and particles free of Pauli blocking. We term this 
pre-payment a “diminished” electronic transition in the case of ]>honon creation (leading 
to Stokes scattering), and an “augmented*' electronic transition in the case of phonon an¬ 


nihilation (leading to anti-Stokes scattering). It is important to note (see also section 6.21 
that it is very difficult to produce an electron-hole triplet (i.e. with a phonon change) 
off-resonance in a mode requiring backscattering. Thus the dispersive modes are more or 
less locked into a given diminishment or enhancement depending on the laser energy. 
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Near resonance, the pseudo momentum q“ is determined by the requirement that the 
energy cost of creating the electron-hole-phonon triplet is just the laser photon energy. 
Typically, most of the energy needed is electronic, on the order of 1.5 or more eV with 
the phonon energy on the order of several 0.1 eV, which is not ignorable. The electronic 
transition energy is given by Ee = E^{q“) — £'^{q“) = huj — £'p/ionon(2q“), where super¬ 
scripts c, u refer to conduction and valence bands. As hui is changed, q" of the phonon 
changes according to the valence and conduction band dispersion, which ideally is a Dirac 
cone structure with light-like linear dispersion of both bands. As q“ changes, the phonon 
energy ©non (2q“) changes according to the well known positive dispersion of about 50 
cin“^ per electron volt for the D band. See figure A. For simplicity, we usually use 
intravalley diagrams even when (in some cases) the process is intervaUey. 

If in a single phonon transition, the phonon is created or destroyed at the time of 
emission, the electron-hole pair formation is not diminished or enlianced. All the initial 
photon's energy goes into the electronic transition, as in DR, and no phonon change is 
yet present. The momenta of electron and hole are created in a matched pair (again, near 
resonance) with pseudomomentum q. The conduction band electron may then elastically 
backscatter q —> — q, in the presence of defects. This allows recombination with the hole 
of momentum q, provided it creates a phonon at — 2q at that precise moment of emission 
(due to coordinate dependence of the transition moment). The emitted light must be 
of the proper Raman shifted frequency to conserve energy. The processes described are 
shown in figure B. Since in this case q did not suffer diminishrnent in absorption, it 
will produce a phonon with slightly larger magnitude -2q (no on the q is present) 
than if produced in absorption with its diminished transition. The Raman shifts are thus 
slightly different depending on whether the phonon was created at the time of emission 
or absorption, since given dispersion the phonon will have a slightly different energy with 
a different q. 

This means that all entries in phonon dispersion diagrams made based on the DR 
model need a small correction, because half of the electron-hole t7ansitions are not the 
one DR supposes, differing instead by a phonon energy. The correction is small, and 
contributes to the broadening (the addition of for example q and q“ peaks) as well as 
slight shift of the center frequency 

The terms in the KHD sum constructively and destructively interfere with one an¬ 
other before the square is taken. The sum includes many different momentum conserving 
electron-hole pairs and electron-hole-phonon triplets leading to the same final states of the 
graphene system. Processes leading to different final states (e.g. a different final phonon 
type or energy) appear in separate, non-interfering terms. 
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4.3 The transition moment for graphene 


Consider the integral involving transition between a valence Bloch orbital at pseiidoino- 
meiitiim q.,,, described in terms of Wannier fiuictions aq^{r — Rq^{^)), and a conduction 
Bloch orbital at pseudomomentuni assuming only nearest neighbor (A with nearest 
B) interactions. The transition moment becomes: 






i=i 




(7) 


is a nearest neighbor vector, i.e. Bj = A + ^j{£). qi-qr Bloch pseudomomenta 
in the valence and conduction band, respectively. Here we have given only the simplest 
form of the transition moment: in fact in our calculations we use density functional theory 
modified Wannier wave fuctions as a function of nuclear* coordinate; see section 

The sum and therefore the transition moment vanishes at the equilibrium position of 
the lattice unless q^, = q^ or — q^. = K, a reciprocal lattice vector, since .F^{qr. A. ^o) 
is the same function of q,. for all A. However, it is not the transition moment at a single 
configuration of the nuclei that is required, but rather the integral over phonon matrix 
elements, i.e. equation Suppose q?. = q^, so we are considering a T point vibration. 
Clearly the exponential is unity. But it is easily seen that A. 7 ^ 0 at 

£ = i.e. ^^(qc. A.^o) is odd at about either G mode. (Moving the A lattice up 
slightly and the B lattice down by the same anrount is a distortion along a G mode T 
point vibration). Thus the transition moment can induce changes by one quantum in the 
G modes, at the F ]>oint. 

A phonon of pseudomomentum k = qv — qc can be induced by lattice distortion 
/ €0 aiKi the transition moment if €k) be- 

comes non-vanishing for q,, — q^ = k or q,, — q^ = k ± K. This happens due to periodic 
undulations in ^^(qr. A, ^k) arising from displacement of the lattice RA(Ck) ^nd RB(4k) 
according to a phonon with wave vector k. However, unless strict conditions are met, such 
phonons, though present, are associated with a conduction band electronic level that is 
Pauli blocked, keeping the phonon silent in the Raman spectrum. Elastic backscattering 
will not help, except for special cases, as described next. If hole doping is present, matters 
can be changed in a fascinating way; see section!^ below. 
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4A Avoiding Pauli blocking 

An important pathway that avoids Pauli blocking begins with qj = —q“: i.e. k = 2q“, 
where the superscript refers to a diminishriieiit of the electronic transition energy by 
the simultaneous creation of a 2q“ phonon (in the case of Stokes scattering); see below. 
The electronic transition will be accordingly diminished by the energy of the k = 2q" 
phonon. The conduction electron suffers a q“ —^ — q“ kick to conserve momentimi, and is 
Pauli blocked. However in this case subsequent defect elastic backscattering can re-align 
it with the hole, allowing recombination. A photon is emitted at a frequency revealing the 
diminished Stokes shift of the 2q” phonon production. Since the diminished q“ depends 
also on the incoming laser energy, the ])honon dispersion will be revealed by changing 
that energy. A 2q“ modulation of the G mode for example will become Raman allowed 
by this mechanism. The modest dispersion of that mode will give rise to a sideband to G. 
which we call G', visible only in the presence of defects. This is exactly what is observed^ 
see figure 


Stokes D phonon produced at photoabsorption 




Stokes D phonon produced at photoemisslon (not the 
same q, E or as when produced in absorption) 




Figure 1: A. Phonon creation in absorption, with diminished electronic energy. Al: 
The electron wave vector q is lowered to q", affecting the jdionon energy through its 
dispersion. The blue transition represents photoabsorptiori with no phonon produced. i/£> 
is the effective frequency of the diminished transition. A2. Defect elastic backscattering 
allows matched electron-hole (no Pauli blocking) Raman Stokes emission. B. Stokes D 
band jdionon creation in emission. Bl. Full absorption of the photon by the electron- 
hole pair. B2. Defect induced elastic 2q backscattering of the conduction band electron, 
allowing Raman Stokes emission and D(2q) phonon production perfectly matched with 
the hole. The energies of the D phonons produced in absorption and emission differ, 
according to q“ versus q and the dispersion of the D band. Anti-Stokes scattering (not 
shown) proceeds very similarly, except augmentation rather than diminishment applies. 


Another important case is q,, = q^ , i.e. k = 0, which was discussed in section 4.3 


This corresponds to creating a F point G phonon that carries no momentiun, but it still 
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Stokes emissiou 






carries energy. Some of the photon’s energy is channeled directly into the G phonon 
energy, diminishing the (in the case of Stokes scattering) electron-hole transition energy, 
including shifting the |q| of the transition nearer to the Dirac cone A' point, as if lower 
energy light had been used: 0.185 eV lower for a 1500 cm“^ phonon. Raman emission is 
active because the electron and hole are born matched in q and really to recombine. Fast 
e-e scattering is the enemy of Raman emission, since if it occurs, the change in electron 
momentum makes Pauli blocking a near certainty. (Stu<lies point to a timescale of a few 
femtoseconds before irreversible relaxation of conduction band electrons by e-e scattering. 
Definitive experimental results [51 [TUI (HI UHl the extremely rapi<l relaxation of 

photo excited electrons due to e-e scattering.) The G band is indeed bright, and benefits 
too from off resonant contributions (imlike D), but the fact that a normally weak overtone, 
2D. is in fact perhaps 10 times brighter than G follows from a fascinating process; see 
sections and below on sliding transitions. 


5 Tight Binding and Density Functional Realization 
of Graphene KHD 

The simplest model Hamiltonian for the single-layer graphene involves only nearest- 
neighbor interactions. However, in the presence of crystal distortions, the hopping strength 
should vary with the pair distance. This is especially important to incorporate here, since 
this will contribute to the dependence of transition moments on geometry changes, needed 
in oiu* KHD theory. 

A two-step ab-initio procediue is used to model realistic hopping parameters. First, 
the density functional theory (DFT) calculations are performed using the Vienna ab ini¬ 
tio Simulation Package (VASP) [HI dS] with the exchange-correlation energy of electrons 
treated within the generalized gradient a]>proximation (GGA) as parametrized by Perdew, 
Burke, and Ernzerhof (PBE)[TSJ. To model the single layer graphene, a slab geometry 
is employed with a 20 A spacing between periodic images to minimize the interaction 
between slabs, a 450 eV cutoff for the plane-wave basis and a reciprocal space grid of size 
19 X 19 X 1 for the 1x1 unit cell. 

Based on the DFT band structure and Bloch waves, the Kohn-Sham Hamiltonian can 
be transformed into a basis of maximally-localized Wannier functions fMLWFlflT] using 
the Wannier90 code. The initial projections for Wannier functions are the atomic 
orbitals and the transformed Hamiltonian is the ab-initio tight-binding Hamiltonian. By 
varying the positions of the basis atoms, the hopping strength t for different pair distances 
r can be extracted and its empirical formula at the linear' order reads: 

t{r) = tiro) + fi{r - ro) (8) 

where tq = 1.42 A, i(7-o) = -2.808 eV and /i = 5.058 eV/A. Figure gives the extended 
tight binding results for the G and G' modes at the top, in two three-panel strips. The 
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Extended Tight Binding for 6,6! and D Modes 



D mode transition moment 
1st derivative at k*q+ K 
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Figui’e 2: (above) The extended tight binding results for the G and G' modes. The valence 
state is at the left. The middle panel depicts the conduction band state in each case, and 
at the top right is shown the derivative of the transition moment with respect to the 
k = 0 G phonon coordinate, and the k = 2q G' (formerly called D'; see below) phonon 
coordinate (bottom right), (below) The first derivative of the transition moment in the 
D sideband mode direction at K + q, K = 27r(^, ^ “ 27r(0, The transition 

moments and their derivatives are integrals over the data in these figures, (left, bottom) 
The first derivative of the transition moment has amplitude of 1.2 in arbitrary units, 
(right, bottom) The second derivative of the transition moment, which has amplitude 
0 . 01 . 

middle panels of each strip depict the conduction band states, and the rightmost panels 
show the derivative of the transition moment with respect to the k = 0 G and the 
k = 2q G' phonon coordinate, shown by the small arrows within the images, (notice 
the imdiilatioris in the atomic displacements). The bottom two panels reveal first and 
second derivatives of the transition moments for the D mode at K + q, where K = 

Consider a transition from A’ + (valence) to K — q~ (conduction), with = 27r * 
(T^ = 27r(0, exactly at the Dirac cone, q gives a small displacement 

from the cone center, a is the carbon-carbon bond length at equilibrium. In this case, 
the electron has a momentum change of and the phonon has momentum 2q“. (For 

compaiison, using the graphene sheet only partly depicted in figm’ej^ the constant part 
of the transition moment has amplitude about 60 in arbitrary units). With displaced 
atoms of amplitude 0.01 A (adjusted according to 2q~ modulation), the first derivative of 
the transition moment has amplitude of 1.2; the second derivative has amplitude 0.01. 

We find that the electronic transition moment has a robust first derivative along any 
choice of the independent and degenerate G mode ]>honon coordinates, and this accoimts 
for their presence in the Raman spectrum. Figure shows the local transition moment 


12 






















along one of those choicas. It is seen to be perfectly repetitive with the unit cell translation 
vector as befits a A: = 0 optical mode. Integration of this 2D local transition moment over 
space yields the transition moment at the given nuclear positions, and repeating this after 
a phonon distortion reveals the phonon coordinate dependence of the transition moment. 

Thus D should be robust (with sufficient elastic backscattering), but simultaneous 
production of 2D seems in doubt. Since 2D is the brightest band, an explanation is 
needed, and is given below in sections and 10 


6 Analysis of Graphene Raman Band Structure with 
KHD 

We now go though a few prominent characteristics of graphene Raman spectrimi dis¬ 
covered and discussed over the years, explained (very directly) by the KHD theory. In 
figure a Raman spectrum obtained by the Hilke group is an average of 60 different 
samples, each with defects, in order to bring out weak bands forbidden in clean, perfect 
graphene crystals. The D band is one such case, while 2D is allowed and bright in pure 
samples. 



Raman shift [cm'M 



Figm’e 3: (left) Averaged defective graphene Raman spectrum on a Si substrate, 
redrawn from reference [THj. (right) Placement of phonon dispersion points found and 
assigned tlirough Raman spectroscopy by Hilke et ai JTH], reproduced with permission. 
D' symbol replaced by G^ appropriate to a sideband of G. 
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6.1 Origin of G band intensity 


The coustaiit part of the electronic trausition moment for arbitrary electronic q is non¬ 
vanishing and responsible for most of the light absorption in graphene. It delivers electrons 
to the conduction band without creating a phonon. A phonon’s creation (or annihilation) 
is the result of a changing electronic transition moment as its coordinate is displaced from 
equilibrium. The more rapid the change in transition moment as a function of phonon 
coordinate, the more likely is the phonon’s creation. 

The G modes have no dispersion: the same q = 0 mode is produced independent of 
laser frequency. The G mode may also be produced in emission. The KHD expression 
has the same transition moment promoting either event. The production of the G by 
either means is a small minority of events in any case, as is any phonon producing a 
successful Raman emission. (It is often stated that typically only about 1 in 10^^ or 10^^ 
incident photons causes a Raman emission). The mechanism based on KHD for G mode 
production in absorption is given in figure 

In section 10.1 we discuss the brightening of the G mode band due to hole doping 
and off resonant scattering. For the G-mode, the phonon production off resonance occurs 
by the same mechanisms familiar in KHD that populate vibrational modes in smaller 
systems off resonance. 


Stokes G phonon 
p' produced at absorption 



Recombination 



Figure 4: Production of a G phonon in absorption, according to KHD extended to 
gi'aphene. 
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6.2 No virtual process for D mode production or others requir¬ 
ing backscattering 

The G mode does not require backscattering to be produced. The electron is ready to fill 
the hole the moment in appears. Not so the D mode. In order to produce a D, a recoil of 
the conduction electron is required. Virtual processes are short lived; the electron must 
fill the hole quickly, leaving insufficient time for elastic backscattering. The off-resonant 
electron and hole cannot remain, since after some time the energy imparted becomes more 
certain; off-resonance, the energy is not sufficient to support this. 

The difficulty is seen in the equation If an off resonance state is created with a 
phonon and qv = —without backscattering this situation arrives Pauli blocked at the 
second transition moment operator, which cannot right the situation without creating 
another, counter-propagating D phonon. Otherwise the numerator vanishes and we may 
say the virtual process does not exist imless the second, counter-propagating D phonon 
is produced. This 2D scenario may indeed take place, over wide zones on the Dirac 
cones, even for the same photon. However among these is a resonant D production and 
subsequent second D emission; the resonant process is expected dominate. Even the 
piurely resonant 2D ]>roduction for fixed phonon wave vector can slide up and down the 
Dirac cone, at a single photon frequency, thanks to the linearity of he dispersion. This is 
our first glimpse of the the transition sliding mechanism presented in section 

6.3 Absence of 2G overtone 

Given the robust strength of the G band, at least a small overtone at 2G would be 
expected, yet the 2G band did not make any appearance in the spectra of Hilke et. 
al. [TS|, where other weak bands were seen for the first time. We have been touting the 
uncommon strength of the overtones in graphene, so this absence must be regarded as one 
of its mysteries and must be explained. There is no group theoretical ban on its existence, 
and like G it is a F point mode, requiring no backscattering. Why is it missing? 

DFT-tight binding calculations of the transition moment show that the second deriva¬ 
tive along the G phonon mode is about 2 orders of magnitude smaller than the first 
derivative. Since the intensity goes as the square of the transition moment, this would 
wash out 2G if it were the only scenario for a 2G. 

However there is another possibility to consider: The robust linear slope in the tran¬ 
sition moment along the G mode could be used twice, once in absorption and the second 
time in emission. Although not a simultaneous production of the phonons, two G mode 
phonons will have been produced, and a Raman band would appear at 2G. The phonons 
would both be pseudomomentum A: = 0. It is easy to see however that the intensity for 
this process must be extremely low: If there is an amplitude of 0.025 (probably much 
too high an estimate, but conservative for this purpose) for producing a single G mode 
in absorption, the amplitude fur two G ]>honons, one in absorption and one in emission. 
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is 0.000625. This corresponds to a probability of two G phonons produced this way some 
1600 times smaller than the probability of a single G phonon production. “One up, one 
down” does not lead to a visible signal. 

Paradoxically, we will see that tliis ''one u]>, one down” mechanism is the key to the 
brightness of the 2D band, and other overtone bands. The 2G mode is unable to increase its 
brightness by one up, one down "transition sliding” (section [^, becaase sliding requires a 
k = 2(] phonon be produced. Sliding greatly amplifies the chance of producing a phonon, 
for example in the case of 2D, or 2G' (formerly called 2D') since many simultaneous 
amplitudes are summed. Transition sliding is a key prinoi]>le made possible by linear 
Dirac cone dispersion (section]^, and a key Raman mechanism revealed here. 

6.4 The G' [former D'] and 2G' [former 2D'] Bands 

Nearby the P point k = 0 phonon, the transition moment can also give rise to k = 
2q“ phonons, giving a momentum kick 2q“ to the conduction band electron. Elastic 
backscattering makes recombination possible, revealing the existence of the sideband to 
G known as D'. As Ferraii and Basko suggested [3], G' would be a better name for D'. 
We adopt this notation, in spite of the checkered history of the G' nomenclature, which 
used to denote 2D not many years ago. 

The healthy 2G' band is the overtone of G' and does not require defect elastic backscat¬ 
tering. Graphene Raman spectra in the literature are often cut off before its ca. 3200 
cm“^ displacement. Its frequency is close to twice that of G'. It seems the 2G' band owes 
its unexpectedly large intensity in the absence of defects to the same transition sliding 
mechanism that benefits 2D; see the next section. 

7 Sliding Phonon Production 

It is a consequence of linear Dirac cone dispersion that the transition making a D phonon 
in absor])tion, shown at the left in figureis allowed over a range of valence to conduction 
electron-hole transitions all for the same incident laser frequency and all producing the 
same phonon. The therefore result in the same final state and contribute to the same 
Rama amplitude, which is squared to give the cross section. 

The brightness of the 2D overtone mode (and of other k = 2q overtone bands like 2G') 
is a consequence of these continuously many simultaneous transition amplitudes, each 
one starting from a different valence band level, all for the same incident laser energy and 
all producing the same D phonon in absorption and one of opposite pseudomomentum 
in emission. Am]>htude for electron-hole and electron-hole-phonon production appears 
simultaneously up and down the Dirac cone, vastly extending the nominal syimnetrical 
A‘ — ^ to K H- q transition (giving a -2q kick to a new phonon). That is, a continuum of 
K — q + Sq to K + q + Sq transitions is available, all generating the same —2q phonon. 
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Reference transition 


Transition upsliding 




Figure 5: Transition sliding on a Dirac cone, (left) The nominal, symmetric '“Double 
D'^ (2D) transition with energies equal above and below the Dirac point. Absorption is 
followed, without any backscattering, by emission along the reverse path, producing a 
second phonon of opposite momentum, (right) By moving the rectangle up, not changing 
any dimension, we find it still fits absorption and emission transitions (it could also slide 
downward) along the Dirac cone. The reference transition is shown as a dashed rectangle. 
The phonon 2q is imchanged, since q is sas is the electronic transition energy and incident 
photon energy. 

This is possible due to the linear dispersion of the Dirac cone. Since the 2q phonons are 
the same, the sliding transition am]>litudes all shai'e the initial state \i) and final state 
I/) and add coherently, before the square is taken. The energy sum over 

electron-hole-phonon triplets r))|Xinvr;(C)) changed by sliding: the same 

photon is absorbed, and always put into making the same —2qv phonon with the balance 
put into the electronic part of the transition; Ev,c,m is highly degenerate. This is true of 
every distinct Ey^cm: resonant or not: a continuum of resonant and off-resonant states 
^))IXm^,c(C)) with noiivaiiishmg transition moments ranges over valence holes and 
conduction electrons (plus 2q phonons) with energy near Ei + fkoj. 

The sliding transitions are already built into the KHD equation one has only to 
examine the possible intermediate states to see that they are present. Their idtimate 
limitations will be a rich topic for future study: trigonal warping and distortion away 
from lineal’ dispersion at energies farther from the Dirac ]>oint will cause slightly different 
energy phonons to be produced; these will add to the intensity non-coherently, only after 
the square. They will broaden the transition but are intrinsically less bright due to 
diminishing constructive interference. Also limiting the sliding is the density of states, 
diminishing to 0 as the band edges ai’e reached. 

We have checked the transition moments for sliding transitions with oiu tight binding 
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model. For example, the transition from an occupied conduction band level at K + q 
to an empty conduction band at K + 3q. The transition moment first derivative is not 
significantly smaller than a transition from an occupied valence band at A" — 5 to an empty 
conduction band at A + 

Sliding adds to the transition amplitude for production of the first 2q phonon, but 
electrons produced by such sliding transitions do not match the holes they leave beliind 
and ai‘e thus Pauli blocked in emission. Elastic backscattering does not rescue these 
phonons from obscurity; {K + q+Sq, backscattered elastically, gives K — q — Sq^ which doas 
not match the original hole at K — q + Sq). Below, we note that hole doping creates empty 
valence bands, that can, if conditions are right, accept emission from such conduction band 
states, giving rise to a spectacular broadband '‘electronic Raman” emission as seen in 
reference [TH]; see section 1(1 (These authors provided a different interpretation involving 
"hotband” emission following electron-phonon scattering). 


7.1 Reversing the path 

The Raman-silent single phonon process becomes Raman active when the conduction 
band electron, produced by a sliding transition along with s. k = 2q phonon, emits (with¬ 
out first backscattering) to the valence band along the reverse path used in absorptiony 
(see figure [^. A second, oppositely propagatmg phonon is released. The electron is 
automatically matched to the hole and recombines, without Pauli blocking of any sort, 
and in the case of D phonons a proper matched 2D phonon pair has been produced. 
The cumulative effect of constructively interfering contributions over a continuous sliding 
range of the terms in the KHD sum (since the numerator of every term in the sum is an 
absolute value squared, before the overall square is taken) causes a large enhancement of 
the 2D Raman band intensity. In graphite whiskers, the intensity of the 2D overtone is 
found to be about 10 times stronger than the single phonon G mode, normally expected 
to be much stronger than an overtone band [20]. As mentioned above, the G mode is not 
amenable to sliding transitions. 

Figure shows the sliding scenario. A large range of sliding AE^' = AE^, i.e. equal 
shifts in valence and conduction bands, are available; all of these are resonant, not virtual, 
transitions. Compared to teh reference, non-sliding symmetric case, in up-sliding the 
valence wave vector has been shortened by 5q, and the conduction wave vector has been 
lengthened by the same amount, so it remains a 2 q“ transition and giving a 2 q“ phonon 
production, just as when AE^' = AE'^ = 0. All the shding transitions are independent 
amplitudes at the same ]>hoton energy simultaneously present, and together they vastly 
enhance the probabihty of producing a pair of D phonons. The density of states for both 
the initial and final electronic states will have a major effect on the propensity to slide 
various amounts. 

The sliding mechanism explains many known facts of Raman scattering in graphene. 
First and foremost, the brightness of the 2D and other overtone bands (mixed transitions 
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can occur as long as the pseudoinorneiitum is the same in absorption and emission) results 
from addition of a continuous range of sliding am]>litudes, before the square is taken to 
give Raman intensity. Second, the fact that the 2D peak strongly decreases with increasing 
doping or disorder is now explained easily: Do]>ing can provide much faster phonon-less 
emission pathways to empty valence states that will diminish production of the emission 
producing another D ]>honon. Elastic scattering of the conduction electron (it does not 
have to be backscattering) by defects will also quench the D emission probability by Pauli 
blocking and thus quench the 2D intensity. 

1.2 The D band 

As we have seen in figiue the electronic transition can slide up and down the Dirac 
cones and still produce a 2y“ D phonon. The conduction band electron aiTives at 5 “ +5q. 
having shifted momentum by 7 ^ 0 and getting a. 2q~ kick from the phonon. It will 
not match its hole at —q~ + Sq, even with elastic backscattering, since it would have 
momentum — 9 ’ — Sq. Thus D phonon production visible in the B.aman signal does not 
benefit from sliding. The D ]>honons are present nonetheless; the electron is hung up in 
the conduction band and must relax by other means than Raman emission. 

7*3 D band Stokes, anti-Stokes anomaly 

The Stokes versus anti-Stokes frequencies in the D and 2D bands are graphene Raman 
anomalies, discussed first in pyrolytic carbon [2110 mid graphite whiskers [22] • There are 
two striking experimental results to explain: (1) A difference between D band Stokes and 
anti-Stokes frequencies. In a small molecule, the Stokes and anti-Stokes bands measure 
the same vibrational state, and there cannot be any difference between the two: they are 
symmetrically spaced across the Rayleigh line, i.e. 0 asymmetry. For the graphene D 
band, the asymmetry is instead about 8 or 9 cm“^. (2) As shown in the next section, 
the 2D band Stokes, anti-Stokes asymmetry is not twice the D band asyimnetry, which 
would be expected because two D phonons are produced, but close to 4 times the D band 
shift, or 34 cm“^. These numbers emerge simply from oin KHD theory, without invoking 
virtual processes. 

The D mode Stokes band is an average of emission and absorption production, with 
emission production imshifted. Production in absorption is shifted down 8 cm“^, by the 
dispersion of D and the diminishment cause by energy conservaton. The average of the 
two is 4 cm“^ closer to the Rayleigh line than the undiniinished emission production 
alone would be (see figure . Similaiiy, the D mode anti-Stokes production in absorption 
also consists of two bands, overall 4 cm“^ higher in energy and fa7iher from the Rayleigh 
band. At 3.5 eV and 1350 cm“^ the Stokes vs. the anti-Stokes D phonons (refiected about 
the Rayleigh line for comparison) will differ by about 8.4 cm“^. This is an anomalous 
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Stokes, anti-Stokes asymmetry of about 8.4 wave numbers, in excellent agreement with 
experiment. Thus the anomaly has a simple explanation in terms of jjhonon production 
in absorption vs. emission. 

7A The 2D band 

The 2D overtone band in pure graphene is the strongest line in the spectrum, even much 
stronger than the fimdamental G band. The Kohn anomaly has been proposed as a 
contributor to the strength of 2D, and indeed it may be, but then G is weaker, also born 
at a Kohn anomaly, and is a fundamental, not a normally weak overtone. 

It is important to note that if two counter-propagating D ]>honons were actually pro¬ 
duced simnltaneously in absorption, there would be a doubling (two phonons) of a double 
diminishment of the electronic energy (since twice the energy is needed from the photon 
to produce both phonons at once). This implies a shift of 32 cm“^ relative to the pre- 
smnably equally important simultaneoiLS emission production of two coiuiter-propagating 
phonons, a transition that is not diminished in energy or q. This would imply that the 2D 
band would either be double or a single peak considerably broader than 32 cm“^ This 
is not consistent with experiments revealing a slightly asymmetric line about 25 cm“^ 
FWHM [231. 

But there is another possibility, just described, that of producing a D phonon in 
absorption and another in a miiTor image emission. There are a continuum of such 
amplitudes which can also slide up and <lown the cone, part of the KHD sum, each 
contributing an imaginary part with the same sign to the total 2D amplitude and each 
producing phonons at the same 2q pseudomomentum (see figiu'e]^. The simultaneous 
2D production in emission or absorption with its 32 cm“^ problem is thus alleviated (not 
to mention it is very weak compared to what sliding produces). The sliding mechanism 
also predicts the experimentally measured Stokes-anti-Stokes anomaly for 2D (see below). 
The density of states for both the initial and final electronic states will have a major effect 
on the propensity to sUde various amounts. 

7.5 2D band Stokes, anti-Stokes anomaly 

The 2D Stokes, anti-Stokes anomaly is simple to explain. In th sliding mechanism, both 
phonons share the same energy and diminishment, or augmentation (in anti-Stokes). Each 
Stokes D (one in absorption, one in emission) suffers an 8.4 cm“^ diminishment, or 17 cm“^ 
total, and each anti-Stokes D enjoys an 8.4, cm“^ augmentation, or 17 cm“^ altogether, 
for an overall 34 cm“^ asymmetry, foiu times greater, not the expected two times greater, 
than the 8.4 cm“^ D band asymmetry. Again, this is in agreement with experiment (see 
figure [^. 
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D and 2D Stokes, anti-Stokes anomalies 

2D scale is compressed by a Solid: Stokes anti-Stokes bands are reflected 

factor of 2 relative to D. Dotted: anti-Stokes about the Rayleigh line. 



Total D Stokes 


Stokes production 
of D in emission 

Stokes production 
of D In absorption 

Total D anti-Stokes 

anti-Stokes destruction 
of D in emission 

anti-Stokes destruction 
of Din absorption 


Figure 6: The prediction of Stokes and aiiti-Stokes Raman spectra of the D and 2D 
band positions and their anomalies, in KHD theory. The horizontal scale of the 2D-band 
spectra (top) spans twice the range of the D scale. Anti-Stokes processes were brought 
across the origin of the Rayleigh line and are shown dotted. See text for details. 

The D and 2D Stokes, anti-Stokes anomalies are thus easily explainable using KHD 
theory, unlike the very elaborate rationale used in the DR model [^. 

7.6 D and 2D bandwidths; 2D intensity with laser frequency 

The D-mode band must be broader than the band centers spaced by 8 or 8.4 cm“^ that 
comprise it. We find it to be ca. 20-25 cm“\ FWHM in the literature, about twice the 
FWHM of G the mode, or 13 cm“^|5]. The G has no dispersion and no issues of Stokes, 
anti-Stokes anomalies. A D line much broader than 13 cm“^ supports the idea that 
the Stokes D mode is really the superposition of two displaced peaks, one coming from 
production in absorption, another in emission. Not only are the Stokes and anti-Stokas 
peaks shifted by 8.4 cm“^; they also must be overlapping when brought to the same side 
of the Rayleigh peak, exactly as seen in experiment; see also figure 

The 2D bandwidth is only somewhat larger than D. approximately 30 cm“^. It is not 
a double peak, at least not until the symmetry is broken (as revealed by tensioning the 
sample in some direction.) It seems likely that 2D earns its width in a different manner 
than D, perhaps a result of the sliding process on slightly nonlinear' or trigonally warped 
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Dirac cones, or off-resonant sliding. 

The decrease in 2D intensity with increase laser energy, and broadening of 2D, leading 
to neai' complete absence at 266 nm in the ultraviolet [25 ^ to explain based on 

the KHD and sliding transition pictm’e. The graphene electronic dispersion is suffering 
significant bending (as distinct from trigonal warping) as the ultraviolet is approached, 
degrading the linearity essential to the sliding 2D mechanism and its interference enhance¬ 
ment. Note that is G mode is robust to the laser energy increase, and indeed it does not 
transition slide. It has been reported to increase intensity at the fourth power if the laser 
frequency [23, which incidentally fits the ojjU)g classic KHD dependence, equation G is 
ready to emit a Raman shifted photon immediately upon photoabsorption, and benefits 
from off-resonance contributions. The first D absorption step in 2D sliding needs to be 
close to on resonance, since the phonon created and the associates electron recoil cannot 
immediately emit, mitigating against a virtual process. 

8 Defect and Laser Frequency Trends 

Several interesting trends develop in the Raman spectra as density of defects or laser 
frequency changes. We begin with the fascinating similarity of sidebands in polyacetylene 
and graphene. High symmetry A; = 0 dispersionless bands can be parents of dispersive 
sidebands carrying momentum 2 g“, coming from the production of a phonon in absorp¬ 
tion, where q~ is slightly less than the higher electronic pseudomomentlun when no 
phonon is produced. The conduction band electron gets kicked to — 5 “ as it generates a 
2q~ phonon. It requires elastic backscattering to appear in the Raman signal. In emission, 
the sideband phonons carry momentum 2qi,. 

The trends in the G' band are quite parallel to sidebands in polyacetylene: fixed A: = 0 
peaks with nearby dispersive sidebands, growing in intensity with increasing soiurces of 
elastic backscattering, and moving in frequency according to the band dispersion and 
phonon q. The G' band (formerly called D') has nothing to do with the D band, and is 
simply the sideband to G, as Hilke [TH] and others have known for some time. 

The growth in sideband intensity with increasing sources of impinity backscattering is 
seen on the left of figure [2S1 • We explained these trends entirely within a KHD context 
applied to the one dimensional polyacetylene crystal with defects j 2 j. 

We reproduce studies of the trends with laser energy in figure]^ left, for both graphene 
and polyacetylene [ 2 H]« hi graphene, the conductance trend is toward increase<l sideband 
strength with decreasing laser energy, as seen in figure [7j In polyacetylene, the trend is 
reversed. One obvious difference is that a propagating electron wave cannot fail to collide 
with any defect in one-dimensional polyacetylene. We have un<lertaken density functional 
calculations on polyacetylene distorted by kinks and other geometrical defects, which show 
that higher energy electronic states backscatter more readily from such defects than do 
lower energy ones, in the conduction band. 
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Defect density and intensity of G and 2D bands 
polyacetylene graphene 
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Figure 7: (top) Response of polyacetylene and graphene Raman spectra to defects, (top 
left) The sideband growth is shown in the 1054 cm“^ polyacetylene band as elastic 
backscattering defects are added. The figure was redrawn and the k = 0 band con¬ 
tribution colored, using Schafer-Siebert et. al [2S]- The sideband is present at 0% added 
defects because of pre-existing defects (including molecular ends). The total intensity 
consists of the A: = 0 band and the overlapping sideband, making a break in slope to the 
right of the k = 0 band peak, (top right) An important graphene Raman study taken 
from Childres et.al^^y with permission, with similar sideband growth (G', sideband to 
G). The dramatic reduction in the 2D intensity with defect density is explained in the 
text, (bottom) Trends in sideband growtl??ind dispersion are shown for polyacetylene 
and graphene (red arrows in the case of polyacetylene, giving band frequency predictions 
based on electron and phonon dispersion [2]). The change with laser energy is opposite 
in the two cases. 


























8.1 The D and 2D bands are sidebands to a forbidden 

k=0 K point vibration 

In a definable sense, the D and 2D bands are dispersive sidebands to a forbidden k=0 
Dirac K point phonon we call living at the vibrational K-point. In analogous cases, 
such as the G' [D'] and 2G' [2D'] bands of graphene (dispersive sidebands to the non- 
dispersive G band), see figure the dispersive sidebands in polyacetylene |2] ai‘e partner 
to visible non-dispersive bands. In those cases the non-dispersive bands belong to F point 
vibrations. The band edge A’ ])oint vibration D^^ is of course not a F point vibration; 
one consequence is that it has a vanishing transition moment and Raman intensity, as is 
required by symmetry. This however does not discjualify the D and 2D bands troni being 
labeled dispersive sidebands to this silent and dispersion-less pai'ent band; it would lie 
typically 8 cni"^ to the left of D and not require backscattering. 

8.2 Defects and 2D intensity 

There is a well known dramatic decrease in the 2D band intensity with increased defect 
density [2^ (figure 0. This is easily understood in terms of the sliding production of 
the two D ])honons. The amplitude for production of the first D phonon in absorption 
is relatively insensitive to defect density. Transitions with no sliding contribute to D 
Raman intensity if they are elastically backseat tered, and indeed the D intensity increases 
with added impurities. Sliding D transitions are Pauli blocked, even if backscattered. 
Sliding transitioned electrons are equally prone to defect elastic backscattering or more 
general scattering, making it extremely unlikely they can produce another D in emission, 
being unable to reverse absorption path. Defect scattering of the conduction electron 
thus quenches the source of D phonon production in emission, and the 2D Raman band 
intensity diminishes with it as defects are added. This is just what is seen in figureThe 
reverse trends of D and 2D intensity with added defects therefore follow from the sliding, 
‘*up and back down along a reversed path*' mechanism for 2D Raman emission. 

8.3 Anomalous spacing of D and 2D 

A prediction can be made about the D, 2D spacing seen in experiments at any frequency. 
This has been discussed within the DR model also [2TI EU]* The frequency of 2D is 
smaller than twice D, by amounts depending on experimental conditions. The ideal 
‘*bai’e,’' unstrained, low temper at m’e, and fairly clean (but dirty enough to see D) graphene 
experiment has not been done to our knowledge. However, quoting Ding et.ai [30], '’The 
results show that the D peak is composed of two peaks, unambiguously revealing that the 
2D peak frequency {uJ 2 d ) exactly twice that of the D peak {ujo )■ This finding is 

confirmed by varying the biaxial strain of the graphene, from which we observe that the 
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shift of CJ 2 D /2 different.*’ 

According to our a]>plioation of KHD to graphene, a 1335 cm“^ D phonon produced 
in absorption is diminished owing to the 50 cm“^/eV phonon dispersion, by 1335/8065 
X 50 = 8.28 cin“^. A D phonon produced in emission is undiminished. The two bands 
will overlap to make a broader feature than either component. Assuming the two bands 
are equally intense, as KHD predicts, there should be a combined band with an average 
4 cm“^ displacement to the left in the Stokes spectriun. The idea that D is composed of 
two bands with an 8 cni"^ splitting was also suggested within the DR model [2T], with 
a very different justification, “'depending on which of the intermediate states is virtual” 
[^. The reason for the two bands is actually much less exotic (absorption vs. emission 
production) and on a firmer foundation in KHD (both real, resonant processes) than the 
virtual processes required in DR. 

The 2D band consists of two separately produced, diminished D phonons. There is 
an 8 cm“^ dirninishment in absorption, and a matched 8 cm“^ diminishment in emission 
according to the sliding scenario, totaling a 16 cm“^ shift. As just discussed the D band 
is displaced by 4 cni“\ thus twice the frequency of D is predicted to be 8 cm“^ shifted 
as opposed to the 16 cm“^ shift of 2D, or a -8 cm“^ difference between 2xD and 2D. 
Review of many published spectra under different conditions shows B{2D) — 2E{D) ~ 
-2 to — 10 cm“h However samples were suspended on different substrates by a variety 
of methods; measiued shifts depend on these conditions. Reference shows that any 
source of stretching or compression can affect the D, 2D distance. The D. 2D shift deserves 
more investigation using suspended, gently pinned graphene. 

9 Mixed Bands and Bandwidth Trends 

Does the sliding mechanism enhancing the 2D band brightness also contribute to the 
strength of other bands? We have already noted that sliding does not help the D band 
gain intensity, since the electron is Pauli blocked even with elastic backscattering (however, 
see the last section for the changed situation when the sample is hole doped). The G band 
transition cannot slide since the P point produces a = 0 ]>honon. The same applies to 
2G, which does not appear in the Raman spectrum regardless of defects. 

The sli<ling mechanism for mixed transitions is a different story than for homogeneous 
ones, especially for bandwidth. The momentum conservation requirements on the pro¬ 
duction of a ])air of phonons requires that they are matche<l in q. They do not need to be 
matched in type; e.g. a Raman band for G' (old D') and (Hilke’s notation) could be 
produced by the sliding mechanism. This fact allows us to explain many of the disparate 
bandwidths seen in the Hilke et.al. data [TH], since differences in dispersion and frequency 
of the two components contributes to the bandwidths, as we presently see. 

Most of the mixed transitions are weak (and certaiidy woidd be invisible without 
transition sliding), as figure shows. The weakness may reflect small transition matrix 
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Overtone and combination band widths 



G’+D7 2G' G’+D3 G'+D" D+D' D+D5 G'+D4 D+G^ 2D+G' 


Figure 8: Calculated (blue) and experimental (from [TS|. r^d) bandwidths for mostly 
mixed combination bands. The calculations were as follows: The wavenumber difference 
for G' (old D') up first and D4 down versus D4 up first and G' down was calculated iLsiiig 
the mode frequencies and dispersion slopes from [IS]. The Raman shift (and phonon 
energies) depend on which transition is first, and this difference is widening the bauds. 
For inodes a and 6, the formula is W = (I*5q — Sb\)ifi^a + f^t) where W is the component 
of the width coming from which mode is created first, Si is the dispersion slope in cm“^ 
per eV, and luji is the energy of the phonon in eV. The theoretical results are compai’ed 
to widths found in the data from [TH] after adding 16 cm“^ for the “intrinsic” width. 
Except for the G' + 2D (assignment and data taken from[^) and G' + D, the agreement 
is good, and it lends confidence to the KHD approach and the one up, one down scenario 
double phonon. See text for discussion of bands labeled by a red star. 

elements or the possibility of destructive interference between terms in the sum in equation 

The strength of homogeneous two phonon transitions, like 2D or 2G', is expected to 
be high since they are produced by mirror image processes with in-phase numerators. 

For mixed transitions, sliding applies but leads to slightly different Raman band fre¬ 
quencies depending on which phonon is produced first, in absorption. The bandwidth will 
reflect this (see figure]^. Using data from Hilke [E] at 288 nm, we arrive at figureThe 
reasons for the bandwidths of the D-h G' and 2D + G' are discussed below. Red stai's ou 
D + D" and G' + D": We have used reference [32] to help imderstand the skewed line at 
about 2450 cm"^ with the nominal assignment D -h D". This study decomposed it into 
two bands, one of wliich is D -h D" with a width of 20 cm“^ and the overlapping higher 
energy LO G' band near the T point, but now near the K point, a less intense band, with 
a FWHM of 29 cm“U 

Some mixed sliding transitions, such as G' + D^, G' + D^, and even some hint of D + 

(Hilke’s notation, except D' -¥ G') do not require defect backscattering and are seen 
weakly in high quality spectra of clean graphene, for example in Childres et. ai, reference 
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Consider a mixed overtone baud involving phonon modes A and B. If A is created in 
absorption, the of the transition will be different than if B is created first, giving an 
electron pseudoinomentmn q^. If A is a higher frequency phonon than B, the electronic 
transition energy diminislnnent is larger if A appears in absorption, and q^ will be smaller. 
The emission B phonon must follow with the opposite pseudomomentum 2q^. This 
allows fur matched electron-hole recombination, but the B phonon is required to adjust 
its energy to arrive at the right q^. This energy correction depends on its momentum 
dispersion. Thus the total phonon energy (and thus the Raman shift) is slightly different 
if A is created in absorption than if B is. This fact contributes to a small uncertainly 
in frequency and contributes to the band width. We calculated the bandwidths for each 
mixed transition and added 16 cm“^ to allow for the intrinsic broadening seen in narrow 
bands, due presumably to phonon decay [311 ESI • Using data mostly from Hilke [E] at 
288 nm. we arrive at figure This gave rise to figure 


10 Evidence of Sliding D Absorption 

Chen and co-workers [E] doped the valence band by as much as 0.8 eV and saw 
abundant continuum emission, in a certain range of depletion and Raman shift. We now 
show this emission, that as reference [TU] points out, integrates to more than 100 times the 
strength of the G baud that it overlaps, can be explained by sliding D phonon transitions 
that are normally Pauli blocked or could have been the first step in making a 2D pair by 
reversing the sliding transition, coiqded with an electronic Raman component. 

With hole doping, D phonons produced in a sliding transition, and the associated 
electronic conduction band level ai‘e normally orphaned by Pauli blocking (unless the 
electron returns by the sliding 2D mechanism along the reverse path). Now the electron 
has a new option: to emit to a valence orbital emptied by doping without creating another 
pho7ion. The original hole remains unfilled, and the valence orbital is filled instead. This 
leads to a continuum of potentially lai'ge electronic Stokes shifts. (Electronic Stokes 
because the valence electron has been promoted, not phonon in the final step.) This 
phonon-less emission channel is far more likely than creating another phonon, and thus 
the ^'feedstock” of the 2D band is depleted, quenching the 2D band. Thus the 2D band 
should fade out in the experiment as the continuum emission appears, just as seen in the 
experiment (see lower dotted white line, figure]^. Before that happens, we note that the 
deeper the hole doping, the less sliding up distance is available before the last populated 
initial valence level is reached. Deeper doping progressively removes sliding 2D transitions, 
weakening the 2D band as 2\Ep\ increases. The consistent, progressive weakening of the 
2D starting as low as 2\Ep\ = 0 is evident in the lower right experimental data. 
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Figure 9: Continuum emission in graphene from normally Pauli blocked sliding transitions 
at some laser energies and hole dopings. At the top, four scenarios for doping levels and 
sliding transitions at 1.58 eV laser energy are shown ^ together with their relation to the 
data of reference [TH]. For a D phonon, the diminished electronic part of the transition is at 
1.41 eV = 1.58 eV (laser energy) - 0.166 eV (phonon energy). Scenario 1 at 2 |F^f| = 1-1 eV 
is Pauli blocked without defect backscattering and therefore invisible in Raman. Scenario 
2 at 2\Bp\ = 1.1 eV is shown sliding down 0.15 eV, thus matching the hole and avoiding 
Pauli blocking. However is produces a minimum 3765 cm“^ Stokes shift (lower dashed 
line). Scenaiio 3 at 2\Ep\ = 1.4 eV is within the bright continuum emission starting at ca. 
1340 cm“^ and going higher, but cut off at higher Stokes shifts by (1) declining density of 
states as the relevant conduction and valence bands approach the Dirac point, and (2) the 
emission factor cjf (see equation [^. Scenario 4 (upper dashed line) needs to slide down 
to reach occupied levels that can be promoted to the conduction band, but this again 
causes large Stokes shifts. Scenario 5 applies to a hypothetical 2.18 eV laser energy (not 
used in the hole doping experiment of reference [H]); it shows that at any available hole 
doping, the D + electronic Raman shift contimuun emission would not appear anywhere 
near the low energy regime seen in the current experiment. The complete fading out of 
the 2D band at the onset of the continuum emission (lower right) as 2\Ep\ is increased is 
explained in the text. The panel at lower right is re-drawn from reference pi?] . 
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10,1 G mode brightening 


In the paper by Wang and his group [IS], a brightening of the G band is noted as hole 
doping is increased. It can be seen as a gradual waxing of G intensity in figui’ej^ even after 
the continuum band is exceeded in the upper left corner of the plot on the lower right. 
The authors attributed this to removal by do]>ing of destructively interfering paths. This 
also happens within the KHD picture. The sum over nonresonant states |n) of energy 
normally extends above and below resonance, wliich causes cancellations in the real part 
of the sum. Doping eliminates part of the sum on the high side of resonance, enhancing 
the real part. The relevant states |n) all contain matched electron-hole G phonon triplets 
and matched electron-hole pairs (if G is to be produced in emission, with slightly different 
energy denominator) even if En is quite non-resonant (see figure p!o| ) 

When present, real (not virtual) processes play a dominant role in KHD, and apart 
from hole or particle doping scenarios real pathways are always available in graphene. 
Virtual processes (such as those present in ordinary off-resonant Raman scattering) do 
not normally play a center stage role, living in the shadows of the real, resonant processes, 
contributing mostly near’ resonance, in accordance with damping factors. 

Of course, off-resonant (pre- or post-resonant) Raman scattering still operates in 
graphene as it does in molecular systems. For hole doped graphene such transitions 
may dominate when the resonant initial valence states are depleted of population. For 
example, starting on the lower Dirac cone below a hole doped, lowered Ep^ a laser may 
be too low in energy to reach resonant levels on the upper cone (see figure 10 middle); 
yet, electron-hole excitation and recombination with no Pauli blocking quickly follow 
upon virtual absorption, on a timescale A# given by the detuning AF* from resonance, 
where A^A£' ~ fi/2. A phonon may thus be created or destroyed, off resonance. Again, 
the coordinate dependence of the transition moment is responsible. Time can become 
too short for electron-])honon scattering or any nuclear motion to develop, even though 
Raman scattering is quite robiLst. Raman intensity must come instead from an instanta¬ 
neous phonon creation/annihilation process, which KHD provides through the transition 
moment coordinate dependence. 


11 Conclusion 

The universally used (except in the conjugated carbon community), 90 year old KHD 
Raman scattering theory has been applied for the first time in graphene Raman scattering 
theory. The results ai’e in excellent agreement with experiment, and new insights (and 
predictions about new hole doping experiments) have resulted. The most important of 
these may be the sliding transitions, explaining the brightness of overtone transitions in 
gi'aphene. 

If KHD is used with constant transition moments (without justification), the resulting 
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No hole doping 



Deep hole doping 


Deep hole doping 




the same G phonon 
is excited in all 
electronic states 


G phonon coordinate-► 


levels connecting to 
blue valence orbital 
off-resonance term 
levels connecting to 
the red valence orbital 
off-resonance term 


would-be resonant 
transition - but 
valence level is empty 


K point 


Figure 10: (left) Non-resonant (and one resonant) electronic processes that can contribute 
to G mode creation. In the case labeled 1, the below resonance electronic transition (of the 
type seen in off-resonance R,aman scattering) transiently promotes the electron from the 
lower blue solid to the upper blue dashed line. There is nothing available (with a G phonon 
produced) well above or below the upper blue dashed state (which is an electron-hole- 
phonon triplet) with a norivanishing matrix element. Process 2 is on-resonance production 
of a G phonon. Process 3 is above resonance with red solid line representing the initial 
valence state and red dashed line indicating the state to which a virtual transition has 
a nonvanishing transition moment. All three cases apply to the same incident photon, 
(middle) With hole doping, some transitions are no longer possible, if the initial valence 
state is empty, (right) The G mode phonon coordinate is displayed as the abscissa, energy 
the ordinate. Filled bands are solid; empty bands are dashed. Tlie density of states is 
indicated qualitatively using the level spacing. The restriction red-to-red, blue-to-blue 
is due to q-conservation (small energy differences due to phonon production not shown). 
A continuum of off-resonant transitions apply for the same laser energy (arrows are the 
same lengths in red, blue, and excluded green examples) and the same phonon. The sum 
over such states is truncated here compared to the no doping case, but the missing terms 
have the opposite sign in their real part, above resonance. Thus hole doping enhances G, 
for reason of removal of destructive interference, in agreement with reference [IH]. Time- 
energy arguments apply: the time spent in the excited state is h times the inverse of its 
energy deficit from resonance. 
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lack of any predicted Raman scattering in graphene might induce one look elsewhere. 
"'Elsewhere” ought to have been to free up the transition moments, but DR instead keeps 
them constant, and furthermore bypasses KHD completely. Instead, ad hoc, one or two 
more orders of perturbatively treated electron-phonon scattering are inserted to finally 
create phonons in the conduction band. The result is a fourth order perturbation theory, 
two orders beyond KHD. It is hard to see why these phonons should not be 100% Raman 
inactive due to Pauli blocking of the radically changed conduction band electronic levels 
after inelastic scattering. Nontheless DR leads to a parallel universe so to speak, if one 
is willing to turn a blind eye to Pauli blocking, that can be used as a model to catalog 
observed phenomena. The underlying physics is entirely different than KHD, in fact it is 
physically incorrect. Since the proven, traditional, lower order (in perturbation theory) 
KHD mechanism works extremely well, the DR narrative must in retrospect be regained 
as an outlier and a model at best. 
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12 Supplementary Material 


12.1 Time domain KHD 


According to the time dependent form of KHD theory the amplitude to scatter 

from phonon state is 


t t' 


= jdt'j dt" - t') Y[t') Gt{t' - f")V(i")G+(i")K‘"') 


(9) 


—oo —oo 

t t' 
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where is the dipole operator polaiization p. We have incorporated a damping fac¬ 
tor r to account for the environmental factors not explicitly included in the Hamilto¬ 
nian. Gq (t' — t") = ^ ^//. G+(f' _ < t" is a retarded 

Boru-Oppenheiiner Green function (one that propagates Born-Oppenheiiner eigenstates 
unchanged except for a phase factor) and V(t') is the light-matter perturbation with 
arbitrary time dependence governed by /3(t"), which we take normally to be exp[ 2 cj/f"] 
corresponding to a cw laser. 

The expression shows clearly that the propagation on the con<luction band Born- 
Oppenheimer potential surface takes place after the transition moment has acted at 
time t" on the initial valence wave function. The transition moment changes the functional 
form of that wave function and the Born-Op]>enheimer Hamiltonian is presented at time f" 
with newly created or destroyed phonons relative to the valence state, before any excited 
state propagation has taken place. This is also clear below after we insert a complete set 
of Born-Oppenheimer eigenstates to resolve the propagator. The excited state propagator 
acts until time t\ when the electron fills the hole, giving the transition moment another 
chance to act. No phonons are created or destroyed during the time evolution in the 
conduction band, according to KHD, nor are they needed to produce the Raman signal. 

**)) = |0(€-r)) |Xm,(€)) the Born-Oppenheimer state before the photon 
interacts. The electronic state ^(€:r) is a function of all electrons at r, depending pai’a- 
metrically on all the phonon coordinates ^ . The dipole moment connecting the initial 
electronic state i and the electron-hole pair state qc.qv- /^qcqv(^)’ ^ function of the 

phonon coordinates defined as 
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where the subscript r reminds us that only the electronic coordinates are integrated. We 
insert a complete set of Born-Oppenheimer eigenstates (they are complete, even if not 
exact eigenstates of the full Hamiltonian) in front of the Born-Oppenheimer propagator 
in equation 

1 = XI l'^q.qv(^:r'))|Xm{^))(X'm(e)K0qc,qv(€;r)|- (H) 

where we have acknowledged by absence of a subscripts on the phonon wave function 
Xm(£) that the jdionon modes do not change upon electron-hole pair formation (in ex- 
tended systems like graphene). Since r))|xni(0) = •£'c,.^m|'^'qc,qv{€:^))IXm(0), 

we have 


t e 






c,u,mvc_ 


oo —oo 




( 12 ) 


Apai't from pre-factors, equation with the insertion equation [TT] can be easily con¬ 
verted (ignoring the second, off-resonant term as usual and gently damping the laser 
field at infinite positive and negative times) to the Raman scattering amplitude for the 
nij —^nif process starting and finishing on electronic state i with incoming light of 
frequency cj/, incoming polarization p, and outgoing a, i.e. equation 

The final state is designated, apart from initial and final polarization, by the initial 
(and final) ground, valence electronic state labeled by i and the final phonon occupations 
labeled by in/. The sum labeled by c, v and ni is over all electron-hole states and phonon 
occupations that connect both initial and final states via the transition dipole D. Here 
reside some surprising and important terms, including the sliding transitions (see main 
text). 

It is illustrative to incorporate the transition moment into the phonon wave frmc- 
tiou ¥^qc.qv;m(€) = {0Xm{0- Phoiioii excitatioiis are included in (the 

''electron-hole-phonon triplets^’) but may be much less common than the ]>ure electron- 
hole pair amplitude. 

Equation can be retinned usefully to a new time domain expression [3S1 • 
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p.(T 

TTij 
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oo 
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dt (13) 


It is important to note that the KHD Raman amplitude is overall 2”^ order, involving only 
perturbation in the matter-radiation interaction. One could go to higher order by adding 
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well known non-adiabatic correction terms to Born-Oppenheimer theory, but we do not 
do that here. Even better, degenerate perturbation theory involving the same correction 
terms might be used to account for Kolm anomalies and possibly other effects. 

The sum labeled by j is over all electron-hole states that connect both initial and final 
states via the two transition dipoles D. Here reside some surprising and important terms, 
including the sliding transitions described in the main text. Equation[T^is useful for many 
things, including miderstanding the effect of a;/ on the Raman amplitude, if something is 
understood about the time dependence of the amplitude qj'm(€)l¥^^ qj espe¬ 

cially for early times. The faster the amplitude grows in time, the more robust it will be 
against u/ lying far from resonance for a given electron-hole state p. This is due to the half 
Fourier transform aspect of the time integral and transient behavior near t = 0. Transition 
moment coordinate dependence permits, at time t = 0, {(p^^ qj q, 

immediate finite amplitude at # = 0. 

Returning to the energy domain, we probe the effect of setting the transition moments 
constant in phonon coordinates: 
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p,c 

m/m, 


(^/) = 
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(i.e. one gets Rayleigh scattering only) 
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This is still 2”^ order, but baiTen of Raman scattering. Something has to be done to create 
phonons. Here we allowed the transition moments to vary, a,s they were always meant 
to, and ai’e often allowed to do in other contexts, these past 90 yeai’s. The DR. model 
left them frozen and went deeper into perturbation theory to generate phonons, which we 
beheve are given a Pauli blocked still-birth with no path of the electron to return to the 
hole. 


12.2 OfF-resonant Raman scattering 

Off-resonance, the effective lifetime in the virtual excited conduction band states is fi 
over twice photon energy gap AE of the laser promotion to electronic resonance 
For very short times well off-resonance, the transition moment is applied twice to the 
initial state as the electron is promoted and then fills the hole; nothing more happens. 
There is not time for nuclear* wave packet motion in cases where equilibrium geometry 
changes in the excited state. In the time domain picture of KHD, wave packet motion 
takes place along Born-Oppenheimer potentials mostly in steepest descent direct ions [57]. 
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electrons are interacting with phonons but most decidedly not inelastically. The displaced 
phonon wave packets are instantly registered as ]>honons as the excited state is reached, 
but the corresponding Raman intensity develops slowly - first order in time - as the wave 
packet develops velocity [S7]. In contrast, the transition moment coordinate dependence 
means instant excited state phonon population and Raman intensity. This time 
dependauce is more robust to off-resonant detiming, and starts to dominate far enough 
off-resonance [33 • Thus the off-resonant D mode contribution suffers even with impurities 
present since there is not sufficient time to backscatter, but G and 2D do well off-resonance 
since no backscattering is required. 


12.3 Bandwidth outliers in figure 

We now discuss the outliers marked by an asterisk seen in figure The biggest deviations 
from the estimates are the D -h G' band, just above 2D, about 50 cm“^ broad as opposed 
to an estimate done with our assumptions of about 25 cm“^, and the 2D -h G', rarely 
reported experimentally at 4280 cm“^ and about 80 cm“^ broad pTT]. Unlike the other 
combination bands, the D + G' band requires impurities and elastic backscattering to 
be seen, allowing it to grow strong (and broad). 2D + G' does not reqiure iinpiuities. 
Presurnabh', different mechanisms are at work in each case, explaining why the bands 
do not fit the assumptions going into the bandwidth estunate and why D -h G' needs 
backscattering. 

Inter valley backscattering is necessary for an electron that produces a D (G') phonon 
upon excitement before it can emit a G' (D) phonon on the way down (see figure Ilf. 


This is because the G phonon has a third of the unit cell of the D phonon, or equivalently, 
the D phonon is at the K point with respect to the G phonon. We fiuther speculate 
that the nature of the scattering (figure allows the D-hG band to be composed of 
both D-hG and G' in the following way: G does not participate in sliding, but it can 
''slant” (the analog of sliding, but for transitions near the vertical). The G becomes a G' 
transition by becoming uonvertical (keeping the electronic energy fixed); this requires the 
creation of low k G' phonons, to keep momentum conserved; a small A:' deviation from 
the electronic K + k is created, which with elastic backscattering of the electron becomes 
a small — fc' deviation from the K — k. The next step, involving emission creating a D 
phonon, varies in energy according to —k' and the D mode dispersion. 

If G Ms produced first, and slides, it is the D emission that must slant on emission, 
causing a range of k values for this dispersive band (although A* = 0 is forbidden by 
symmetry). Between these two possibilities, there is easily a sufficient energy range of 
phonons thus produced by slanting, accomiting for both the enhancement of the “D-hG” 
band intensity with backscattering defects (due to slanting) and the breadth of the band. 


The same slanting transitions can occm* when producing lone G phonons in absorption, 
so why is the G band not corres])ondingly broadened? The reason is that the resulting 
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D + G' slanting transitions 


2D 




Figui’e 11: Genesis of the 2950 cm“^ band labeled, variously D+G (sometimes) or 
D+G'(usually). It requires impurities to become visible. If D is produced by a sliding 
intervalley transition, then after elastic backscattering the G' transition slants to reach 
the hole, generating a (generally) low q G' phonon to preserve momentum. The slanting 
transitions have somewhat larger Raman shift than pure G; the totality of such transi¬ 
tions explains the width of the 2950 cni“^ band. The required backscattering explains 
why impimities are required to make the Dh-G' combination overtone appear, unlike the 
other combination overtones, which can be faintly seen in a relatively pure sample at the 
top right of figure 6 of the main text, Ne = 0, but ai’e suppressed as impurities are added. 
(The 2D-I- G' does not require iinpiuities to be robust [3T].) 


low k G' ]>honons are Pauh blocked, and elastic backscattering only blocks them fiuther. 
The G bandwidth may however reflect the intrinsic '^Pauli blocking tolerance’^ for very 
small k'. 

The only tliree phonon band we discuss in this paper is G + 2D (or G' + 2D), which 
is rai’ely reported experimentally at ca. 4270 cm“^ Raman shift [21]. It is broad, with 
roughly an 80 cm“^ line width, and does not require impiuity backscattcjring. The reason 
for its existence and its linewidth have a plausible scenario from our KHD based approach, 
including transition sliding. 

The band may be produced by first creating a G' phonon in parallel with a D in 
absorption, with a 2q momentum kick to the electron, followed by a -2q emission along 
the reverse D path creating a second D phonon. The new twist here is that the first 
2q kick can be shared in any proportion as 2q = 2qQ' + 2qp). Again, even without a 
detailed calculation of intensity distribution and fall-oif, it is clear that there is more than 
80 cm“^ energy difference available depending on the ratio of G' to D in the first step. 
The emission is through a normal D at momentum 2ci. The 2q = 2qQ' + 2qj^ process 
can also happen in emission. 
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12,4 The “molecular approach” 


All important prior and non-DR perspective on Raman scattering in graphene has been 
termed the ''molecular approach” [SSlEHl* A molecular* polai*izability context of the type 
familiar* from off-resonance Raman scattering was used. Electron-phonon '‘scattering” 
plays no role. The off-resonance Placzek polarizability derivatives for computing Raman 
intensities, like the experimental Raman bands, are often quite similar to on-resonance 
spectra for* the conjugated hydrocar*bons, except for* overall intensity. A finite Placzek 
polarizability derivative requires non-constant transition moments. However, an approach 
restricted to off-resonance polarizability cannot be regarded as a complete theory for 
resonance Raman scattering in graphene, but it is a step in the right direction away from 
DR. 

Consistent with its off-resonant char*acter, the polarizability approach is a near in¬ 
stantaneous picture, leaving little to no time for electron-phonon scattering in the excited 
state. The D mode requires elastic backscattering (and some time) to become visible in 
the Raman spectrum, and indeed the polarizability pictine is most successful it seems 
with the G and 2D modes, which require no backscattering. The molecular polarizabihty 
picture is a big step in the right direction, quite distinct from DR methods. 

Another interesting and instructive contribution to the molecular approach is found 
in Tommasini et.al. [20]. attempting a more general electronic resonance formulation. 
It was based on a KHD fomidation but still resorted to the Condon approximation, i.e. 
constant transition moments. 
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